Josephson junction with magnetic-field tunable ground state 
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We consider an asymmetric 0-n Josephson junction consisting of and tt regions of different 
lengths Lo and L,r- As predicted earlier this system can be described by an effective sine-Gordon 
equation for the spatially averaged phase ip so that the effective current-phase relation of this system 
includes a negative second harmonic cx sin(2i/'). If its amplitude is large enough, the ground state 
of the junction is doubly degenerate ip — ±</?, where ip depends on the amplitudes of the first and 
second harmonics. We study the behavior of such a junction in an applied magnetic field H and 
demonstrate that H induces an additional term oc H cosip in the effective current-phase relation. 
This results in a non-trivial ground state tunable by magnetic field. The dependence of the critical 
current on H allows for revealing the ground state experimentally. 



Josephson junctions (JJs) with a phase shift of tt 
in the ground state [ij attracted a lot of interest in 
the recent years These JJs can be used as on- 

chip phase batteries for biasing various classical Q and 
quantum circuits. This allows for removing exter- 
nal bias lines and reducing decoherence. Currently, it 
is possible to fabricate simultaneously both and tt 
JJs using various technologies such as supcrconductor- 
fcrromagnct heterostructures[l3| or JJs based on d-wave 
superconductors [ll - 13 1 . 

It would be remarkable to have a JJ (a phase bat- 
tery) providing an arbitrary phase shift (p, rather than 
just or TT. Long arrays of O-tt-O-tt-. . . JJs with short 
segments were sugges ted as systems, where the tp J J 
can be realized |1 44 ld |. Each of the segments are as- 
sumed to have a standard current-phase relation (CPR) 
js = ±jo sin 0. The phase cf) in these systems can be writ- 
ten as a sum of two terms 4>{x) = ip + ^(a;) sin?/;, where 
■0 is a constant spatially averaged phase and \£,{x)\ ^ 1 
is alternating on a scale of the order of the segment's 
length. The effective CPR then reads [HOl 



js = ji sin(V') + j2 sin(20). 



(1) 



where ji = {jc{x)) is the spatially averaged critical cur- 
rent density and j2 < is the amplitude of the effective 
second harmonic. We note that special stochastic dis- 
tributions of facets may also lead to j2 > [l7|. The 



value of j2 depends on parameters of the junctions. For 
j2 < — ji/2 the ground state of the system is doubly 
degenerate with tJj = ±ip, where 



if 



-Ji/(2j2)] 



(2) 



This (fi Jj[l6[ is the generalization of a tt JJ and can 
provide an arbitrary phase bias < (p < tt. These JJs 
have unusual physical properties [18|, e.g., non quantized 



Josephson vorti^s that were predicted[lj] and observed 
experimentally P^j . 

The simplest system, which is suitable for the realiza- 
tion of a "short" JJ, is an asymmetric (i.e., Lq ^ L^^) 



O-TT J J with relatively short segments of length L^, L^^ < 
Aj, where Xj oc I/a/to is the Josephson length. Such a 
system can be considered as one period of an infinitely 
long O-TT-O-TT-. . . chain. The best method to reveal the tp 
JJ experimentally is to include it in a superconducting 
loop (thus forming a SQUID) and to measure the sponta- 
neously generated flux. Another option is to study these 
JJs in a magnetic field H. It is worth mentioning that 
there is an apparent contradiction here. On one hand, the 
usual linear phase ansatz (t>{x) = Hx + 0o in asymmetric 
0-7T J J results in a critical current Ic{H) with a cusp- 
like minimum at H = 0, see Eq. (11) in Ref. [Uj as well 
as Refs. 3, 2^ This Ic{H) should be mirror symmetric 
with respect to the Ic axis, i.e., Ic{—H) ~ Ic{+H). On 
the other hand, the presence of the second harmonic in 
the CPR H]) and the linear phase ansatz for leads to a 
Fraunhofer-like Ic{H) with the maximum at H = o[3l. 
Moreover, when \j2\ > 0.5ji there are two critical cur- 
rents Ic±{H), corresponding to depinning of the phase 
from different potential wells, but both of the Ic+ (H) 
dependences have the main maximum at H = 0.[l8| 

In this paper we study the asymmetric O-tt JJ in a 
magnetic field in detail. We demonstrate that the ampli- 
tude of the second harmonic j2 in Eq. ([IJ depends only 
on the JJ parameters, while the magnetic field H induces 
an additional term oc H cos ip. Further, we predict the 
dependence of Ic{H) at small magnetic fields. 

Consider an asymmetric JJ shown in Fig. [T] The de- 
pendence jc{x) is given by 

jcix) = -jo, if - Lj, < X < 0; (3a) 
jc{x) = +jo, if < a; < Lo (3b) 

The behavior of the phase (j){x) is defined by the static 
sine- Gordon equation 



(f>"ix) - jcix) sm[(j>{x)] = -7, 



(4) 



where the coordinate x and lengths Lo.tt are normalized 
in the usual way to Xj and 7 = j/jo is the normalized 
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FIG. 1: (Color online) Geometry of the O-tt JJ. 



bias current density. jc{x) below is also normalized to jo- 
Next, we write 14 1 



jc{x) = (jc) [l + 9{x)], 



(5) 



where 



{jc 



1 



jcix) dx = jo 



Lq — 



(6) 



is the average value of the critical current density and 

-2Lo 

l{x) 



Lq — L-ij 
+2L^ 



for a; < 0; 



(7) 



Lq — L-n 



for a- > 0, 



is the deviation from the average value {{g{x)) — 0). 
Thus, we rewrite Eq. ^ in terms oi g{x) and arrive to 



0"-Oc) [l+g(x)]sin</) = -7. 



(8) 



In what follows we treat short JJs {Lq,L.,^ ^1)- There- 
fore, the solution (f){x) can be sought in the form 14 1 



{x) = ip + ^(x) sin?/;. 



(9) 



where ij) \s a constant j2a| and £,{x) describes small varia- 
tions of the phase around V', i.e., \S,{x)\ <^ 1, {S,{x)) = 0. 
Substituting the ansatz © into Eq. ([5]) and expanding 
to the first order in £,{x) we get 

^"sinip- {jc) [l+g(a:)][l-|-^(.T)cos?/']sin?/' = -7. (10) 



Equation (|T0| has two types of terms: the constant ones 
and varying ones. Note, that the term g{x)S^{x) has both 
the constant (average) part {g{x)^{x)) and the devia- 
tion from the average g{x)^{x) — {g{x)(,{x)) . Thus, from 
Eq. (Hn]) the relation for the constant terms reads 

7 = (jc) [sintp+ {g{x)£_{x)) sin ■0 cos ■0] . (11) 
The equation for (,{x), recalling definition ([S]), is 

e - 3c{x) cosV'e(x) = {jc) [g{x) - {gO cos V] , (12) 



It turns out that both terms cx cos ip have an extremely 
weak influence on the results. Therefore, for the sake of 
simplicity, we omit them right away to arrive to 



r - (jc> 9{x). 



(13) 



Solutions of this equation include four constants to be 
determined from the matching conditions at x = 

e.(o) = eo(o); c;(o) = e^(o), (i4) 

and boundary conditions at x = —L^^ and x = Lq 

d- Lt,) sin ip ^ h; (,'„{Lq) smi/; ^ h, (15) 

where h = 2H / Hd is the normalized applied magnetic 
field, Hci = ^qI {ttXjK) is the penetration field and A « 
2 A/, is the effective magnetic thickness of the J J. Thus, 
we arrive to the expression for £,{x) and can calculate 



(.90 ^ro+F;,^, 



where the coefficients 
To ^ 



4 LlLl 



LqL^ 



3 Ll — 



Lq — Lt^ 

Thus, Eq. (dH) gives the effective CPR of the J J 



7 = (jc) 



sin + Thh cos + s\i\{2'ip) 



(16) 



(17) 



(18) 



The result is remarkable — as in earlier works 14 - lij 
one gets the second harmonic sin(2'0) with the negative 
amplitude of ro/2. In addition, the magnetic field results 
in a ThhcosTp term which additionally modifies the CPR 
and is tunable by magnetic field. Fig. [2ja) shows the 
effective CPRs 7(1/') of a O-tt JJ for several different values 
of magnetic field h. At /i = we might have a doubly 
degenerate ground state ip = ±ip which, upon application 
of h, transforms into a single one. 

The Josephson energy corresponding to CPR is 



U{^) = {jc) 



1 — COS ip + Thh sin + — sin ijj 



(19) 



The plots of U{ip) for the same set of h values is shown 
in Fig. [DJb). Following the evolution of the curves at 
different h one can see how the doubly degenerate ground 
state transforms into a single one. Moreover, one can 
see that at non-zero h the potential energy U {ip) lacks 
reflection symmetry and therefore can be used to build 
Josephson phase ratchets 21, 2^ . 

The effective CPR (|18p at /i = allows to calculate 
the domain of existence of nontrivial solutions, which is 



defined as |Fo| > 1 [18[. This results in 

L-K ^ Lq\\ -T-r^ -\ Lq > Lt, 



4L2-f3' 



4L2+3- 



(20) 
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FIG. 2: (Color online) (a) effective CPR 7(V') and (b) effective 
Josephson energy (7(-i/)) of a O-tt JJ with Lo = 1, L,r = 0.9 for 
several different values of magnetic field h given next to each 
curve. The arrow in (b) points to a local energy minimum, 
which appears at certain values of bias current. 



Since we have made certain approximations, namely used 
Eq. ((T3|) instead of Eq. (fT2|) . it is worth to compare the 
domain ([20|) with the exact resuhji^: which reads 



L-n > arctan(tanh(Lo)); 
Lq > arctan(tanh(L7r)), 



(21a) 
(21b) 



for our case j° = j^. Both boundaries, the exact one 
given by Eq. (|2ip and the approximate one given by 
Eq. (PD|) are shown in Fig. [31 It is seen that our ap- 
proximation works extremely well not only for small Lq 
and Ltt, but also in the limits Loj^-n- — ^ oo, where it 
deviates from the asymptotic value by only ^ 10 %. 

In the J J with the CPR ([T51) and energy one may 
have several stable static solutions and several critical 
(depinning) currents corresponding to the escape of the 
phase from the relevant energy minimum. To find the 
critical currents for a given magnetic field, we look for 
extrema of the CPR (|T8|) with respect to ip. We arrive to 



(22) 



cosV' - T/j/i sin -0 - ro(2cos^ ip - 1) ^ 0. 



This problem can be reduced to a solution of a fourth 
order polynomial, i.e., all solutions of Eq. ((22|) can be 
found numerically. As a result we obtain up to four rele- 
vant roots and up to four corresponding critical currents. 
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FIG. 3: The domain of existence of a state {h = 0) 
on the LqjLtt plane. The exact boundary (|2ip (continu- 
ous lines) and our approximation (|20|l (dashed lines) agree 
very well. The domain has regions of different colors from 
green (approximation |^(a;) sin7/;| <C 1 works well) to red 
(approximation is invalid). The boundaries correspond to 
|CmaxSin(^| — 0.1, 0.2, 0.4, 0.8. The dots indicate parame- 
ters used for discussion of Fig. 2] Vertical and horizontal 
dash-dotted lines show asymptotic behavior of exact bound- 
ary 



Several examples of 7c (/i) dependences are shown in 
Fig. m One can see that for the parameter set (a), which 
corresponds to a state deep in the (p region at h = 0, see 
Fig- El one observes a characteristic rotated diamond-like 
shape with four critical currents in total for \h\ < 0.6. 
The two branches that meet at ft, « +0.6 (—0.6) corre- 
spond to the escape from the left (right) potential well 
in Fig. mb). The upper (lower) branch corresponds to 
the escape in the right (left) direction. At « 0.6 the 
branches meet indicating the disappearance of the cor- 
responding local energy minimum. The two critical cur- 
rents 7c± (in each direction) observed at h ~ coinside 
with 7± obtained in ref. [isl One may be able to ob- 
serve the lower one experimentally in a system with low 
enough damping using a special sweep sequence [18|. For 
some values of h the smaller minimum can even be absent 
at 7 = [see the arrow Fig. HJb)] and appear at larger 
I7I [the vertical dotted line corresponding to h = 0.5 
in Fig. Hl^a)]. The small energy minimum appears for 
0.03 < 7 < 0.06, which also be seen in Fig. Mjo). An- 
other distinct feature, which might be measurable in ex- 
periment, is the shift of the main minimum to h = ±0.11 
in Fig. mja). This shift is the evidence of a non-trivial 
phase state in the junction. Finally, at large \h\ the up- 
per branches of jdh) approach the asymptotic lines 



(23) 



which are also shown in Fig. [H We remind that our 
analysis is valid only for small values oi \h\ and does not 
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FIG. 4: (Color online) The critical current 7c(/i) for O-vr JJs with different segment lengths. Dashed lines show the asymptotic 
behavior (|23p for large \h\. Gray (pink) solid lines show the results of direct numerical simulations. 



reproduce the global "fdh) features, presented elsewhere. 

When the O-tt JJ gets more asymmetric, see Fig. IH^b) 
and (c), the lower jc{h) diamond-shaped domain be- 
comes thinner and finally collapses. The point of collapse 
corresponds to the crossing of the (p-domain boundary in 
Fig. [3] Looking at 7c(^) one can see that the main cusp- 
like minimum shifts with asymmetry which can be used 
to extract the asymmetry from experimental data. 

We also have performed numerical simulations of the 
7c(/i) dependence using StkJJ software [i^. The above 
analytical results coincide with the results of direct nu- 
merical simulation, shown by pink (gray) in Fig. |4l with 
good accuracy. This is especially so at ^ 1 when the 
approximation |^(a;)sin-0| <C 1, used to derive ([TOl) . is 
valid. 

To summarize, we have shown that the effective CPR 
(|18p includes the term cx HcostJj, which allows to tunc 
the ground state. Thus, our system is a hybrid between 



ip pj[ and (fio [25| JJs introduced earlier. The corre- 
sponding Josephson energy profile can be made asym- 
metric allowing to build Josephson phase ratchets. We 
have clarified how the Ic{H) dependence of asymmetric 
O-TT JJ looks like for low magnetic field. It does have a 
minimum, but the minimum is shifted to some zt/imin for 
positive and negative bias, respectively, so that the Ic{H) 
curve is point symmetric. 
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